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ABSTRACT 

The study focuses on the modeling of the temperature profile during the 
fermentation of beer and the selection of the modelled temperature to 
simulate the growth of a microorganism, the consumption of proteins and 
the formation of aromatic compounds (Ketone and esters). The objective of 
the study was to determine how to select the best temperature for beer 
fermentation and how a portion of the biochemical reaction occurs with the 
controlled (selected) temperature. Finite element modelling has been used 
for heat transfer modelling and COMSOL Multiphysics version 5.3 has been 
used for implementation. Version 17 of MATLAB was used to simulate 
biochemical changes with the chosen temperature. The simulated results 
showed that at high coolant flow, a low-temperature profile was recorded 
over the fermentation time. As such, the observed temperatures were 
1.2m3/hr, 1.3m3/hr and 1.6m3/h, 20 oC, 18 oC and 12.5 oC, respectively. 
The modelled vorticity results also indicated that at a flow rate of 
1.2m3/hr, there was a consistent flow of liquid around the agitation center 
relative to other coolant flows. Isoleucine was exhausted after 13hr at 
12.5°C, 80hr at 18°C and 16hr at 20°C from the start of fermentation. The 
simulated results also indicated that ethyl acetate had reached a hold-back 
value of 0.114mol/m3 at 70hr at 12.50C, 30hr at 180C, and 22hr at 200C. 
However, isoamyl acetate retained a retention value of 0.0105 mol/m3 
until the initial concentration of sugar and amino acids was exhausted 
(throughout fermentation) at all selected temperatures. Valine decreased to 
nearly 195hr at 12.5°C, 120hr at 18°C and 85hr at 20°C. The simulated 
nutrient results were again shown to be zero in 210hr at 12.5°C, 110hr at 
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18°C and 90hr at 20°C of luicine consumption. 
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1. INTRODUCTION 

Fermentation produces energy from the oxidation of 
organic compounds. In the case of beer, fermentation is an 
important process unit operation. The initial chemical 
reaction pathway is the conversion of sugar into ethanol 
and carbon dioxide molecule, coupled with exothermic 
reaction and heat and biomass growth (MacDonald et al., 
2008). Concurrently, a number of organic compounds 
(many of which contribute to or compromise beer flavor) 
are formed at low concentrations, due to the multitude of 
side reactions. 


Beer fermentation is one of the large-volume industrial 
processes that are essentially controlled manually. There 
are several reasons that excluded automation from this 
process so far. One is that the essential process state 
variables are difficult to measure on-line. Another is that 
the fermentation process appears to be a complex system, 
which had not been modeled precisely enough so far. On 
the other hand, however, the brew masters make excellent 
beers with their current technology and any automatic 
process supervision and control system must compete 
with these highly skilled human controllers (Sablayrolles 
JM and Ball CB, 2004). Automatic systems would only be of 


interest to breweries if they would help the brew masters 
to produce the same quality of beer in a much less 
expensive way. There are some obvious entry points for 
further auto-mat ion in beer fermentation. In beer 
breweries, many tasks presently being performed by men 
as well as women, such as manual cleaning, filling, 
pitching, measuring, temperature’ control, yeast 
harvesting, pumping into lager vessels, etc., could be 
automated (Gaurdia, 2001). 


The fermentation of alcohol is the most time-consuming 
step in brewing. It can take up to two weeks to complete. 
Fermentation progression is sensitive to yeast pitching 
rate, dissolved oxygen content, batch pressure and 
temperature. The system temperature strongly affects 
yeast growth and metabolic rate: as long as yeast cells are 
kept below 30°C and not damaged, higher temperatures 
accelerate fermentation (Skye&Ling, 2012). Beyond this 
temperature, ethanol and volatile flavor compound loss 
rates are too severe, coupled with increased production of 
undesirable substances and bacterial growth. Brewers 
control ferment or temperature within a narrow range 
during batch progression, to accelerate the fermentation 
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while also ensuring that yeast is not deactivated due to 
denaturation and that no undesirable flavor compounds 
are produced. 


2. Model governing equations and description 
Diffusion and reaction of temperature flow in the beer 
fermentation tank is integrated to an ordinary differential 
equation as written below. It can be solved by finite 
element method. The finite element method is a numerical 
approach by which general differential equations can be 
solved in an approximate manner. The differential 
equation or equations, which describe the physical 
problem considered, are assumed to hold over a certain 
region. This region may be one-, two- or _ three- 
dimensional. It is a characteristic feature of the finite 
element method that instead of seeking approximations 
that hold directly over the entire region, the region is 
divided into smaller parts, so-called finite elements, and 
approximation are then carried out each element. The 
collection of all elements is called a finite element mesh. 
When the type of approximation which is to be applied 
over each element has been selected, the corresponding 
behavior of each element can then be determined. This can 
be performed because the approximation made over each 
element is fairly simple. 


During beer fermentation conduction-convection mode of 

heat transfer occurs from between beer film and bulk 

beer, between bulk beer and the coolant. Equations which 

govern this heat transfer are described as follows. In beer 

fermentations the governing equation for flow of heat 

between medium and the growing yeast is derived from 

the energy balance as follows. The basic assumption of the 

model governing equation are; 

> Cooling jacket is well insulated. No heat loss to the 
surrounding. 

> Specific heat capacity of the wort and the coolant 
remains unchanged over entire fermentation time. 

> Density of wort and the coolant remains unchanged 
over entire fermentation time 

> Assume wort is homogeneous liquid. 

> Assume the flow is laminar 


Energy balance for fermentation tank 
Rate of energy accumulated = rate of energy in- rate 
energy of out + rate of energy generated- rate of energy 


lost. 


< = 0-—0+ Egen — Q,no rate of energy in and out. 


Q =U*A*ATim 


Hence 
d(mC,T)/dt = -X(AHpc*1)*V-U*A*ATim, WhereATim= T-T;/ 
InT/T; 


m = pV 

d(p*V*Cp*T)/ dt = -X*AHge*p1*V-U*A*ATim 

dT /dt = -X*AHgc* pi / p*Cp - U*A* ATim/ p*Cp*V (1) 

Energy balance for cooling jacket 

Following the same procedure from which is energy 
balance on cooling jacket; Rateofenergyacumulated = 


Rateofenergyin — Rateofenergyout + 
Rateofenergygenerated — Rateofenergylost 


dE 
om Hin — Four +9 + Q 


d(m;*C,.*dT;) /dt = Oc*Fe*Cpc* Te — Oc*Fe*Cpc* Tj +U*A* ATi 
mM = Oc*V; 


Hence 
dT/dt = Oc Fe*Cpc* Te / Oc*Vj}*Cpc a Oc Fe*Cpc*T;/ Oc*Vj*Coc 7 
U*A* ATim 


The governing equation for heat flow between film beer 
and the wall is the conduction-convention heat transfer 
which is; 


OT OT 0*T | 0*T 
ply (us + 0S) = k = 


eh oy) in vector form as 


follows; 
dz*p*Cy—+dz*p*C,*uVT + V.q = dzQ + qo + 
dzQ, + dzQ,(3) 


For laminar flow the incompressible Naiver stokes 
equation and continuity equation is used. 


Et — + 4 ,xmomentumequation 


V 


ou Te, SB) —V)er 4. ,, (28 4 Ov 7 
p(w tv) = tHe (Sst a) + PGB 
Tref ), ynomentumequation 


In vector form the two-equation combined as follows. 


mn + p(u. V)u = V[-pI + u(Vu + (Vu)")] + F(A) 


du. a _ 4 . 
‘i + rome contuinityequation, but in vector form 
pV.(u) = 0 (5) 


For transport of diluted species; 


<= + V.(—DiVCi) + u.VCi = Ri(6) 


Ni = —DiVCi+ uCi 


2.1. Model geometry 

Most of the beer processing industry uses cylindroconical 
Shape of fermentation tanks. Therefore, the selected 
geometry was 2D axial symmetry of cylindro - conical 
Shape with mixing direction and flow direction from upper 
to lower as shown in figure 2.1. 





Figure 2.1 Model geometry for finite element 
modeling (Yusuf and Muray Moo, 2003). 
2.2. Nutrient consumption kinetics model 
(Warnasooriya, 2011) 
Formation of some of flavor compound such as fusel 
alcohols have related to consumption of certain amino 
acids. These amino acids are utilized by the cell primary as 
structural components in the construction of structural 
and functional proteins. The rate of amino acids 
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assimilation by yeast is therefore modeled as negatively 
proportional to the yeast growth rate and this is limited by 
availability of amino acids in the media. 


Luicine consumption rate 
dL ax L 
* D (7) 


— = —Y);y* — * 
dt LX" dt) K,4+L 





Where Y\x = Yield coefficient moles luicine need per mole 
biomass growth, K, = Michaels constant for luicine 
(mol/m/3), L = Luicine (mol/m4), D = First order time delay 
(Michiko et al.). 


Isoluicine consumption kinetics rate 


dl dX 1 
— = —Vy*x — *« — *« D(8 
dt a aes, (8) 


Where Yix = Yield coefficient of isoluicine needed per mole 


biomass growth, K; = Michaels constant for isoluicine 
(mol/m/%), D = First order time delay (Michiko et al.). 


Valine consumption kinetics rate 
dV ax V «D (9) 


— = ~—Yvy* — * 
dt VA" dt Ky+V 





Where Yvyx = Yield coefficient of Valine needed per mole 
biomass growth, Ky = Michaels constant for valine 
(mol/m3), V = valine (mol/m/). 


In all cases first order time delay is represented by: 
D=1-e tr 


Yeast growth kinetics 
dx 


7 Hx * X (10) 
Where 
Ux=(Yxo*H1t Yxu*U2t Yxn*U3 ) * 


of yeast growth. 


Kx : epi 
——_——— , is the specific rate 
Kx+(x—-x0)2 


__** ____ ig the indication of cell growth inhibition of 
Kx+(x—-x0)*2 


unsaturated fat. 


3. RESULT AND DISCUSION. 


Where Yxq = Yield coefficient mole X per mole G, Yxm = 
Yield coefficient mole X per mole M, Yxn = Yield coefficient 
mole X per mole N, Kx = Empirical yeast growth inhibition 
constant (mol/m?)2, X, = initial yeast concentration 
(mol/m+?) 


Esters model(Warnasooriya, 2011) 

The synthesis of aroma active esters by yeast is of great 
importance because they represent the largest group of 
flavor active compound in beer (ethyl acetate, isoamyl 
acetate, ethyl caproate and etc.). 


Ethyl acetate formation kinetics 


— =Ypa*(iitpotps)* X (11) 


Where Yea = Ethyl acetate yield per mole sugar fermented 


Where pl = a 

KG+G 
1), Kg = Consumption rate constant of glucose (Mol/m+), G 
i *—_Where um = Specific 
uptake of maltose (h“‘), M = maltose (Mol/m?), Km = 
Consumption rate constant of maltose (Mol/m/%), K’¢ = 
Consumption rate constant of glucose (Mol/m3), G' = 
HNN KG . KM 
KN+N KG+G KM+M 
Specific uptake of malt trios (h!), N = Concentration of 
malt trios (Mol/m3), K'y = Consumption rate constant of 


maltose (Mol/m+). 





Where uc = Specific uptake of glucose (h- 


= glucose (Mol/m?),u2 = 











glucose (Mol/m?).u3 = Where Un = 


Ethyl caproate formation kinetics. 
dE 


tayo * 

pot Yec*ux*X(12) 

Where Yec -Yield coefficient for ethyl caproate per mole 
biomass growth, ux = specific rate of yeast growth (h*!). 


Isoamyl acetate formation kinetics. 


<= Yiac*piac*X(13) 





Where Yiac = Yield coefficient for mole Isoamyl acetate 
produced per mole Isoamyl alcohol formed, IAc = Isoamy] 
acetate concentration (Mol/m+). 


3.1. Effect of coolant flow rate on fermentation temperature and flow condition. 

The temperature profile, flow conditions and vorticity of flow have been discussed in this section at a different coolant 
flow rate. From the modeled results, the temperature profile of beer fermentation was observing with fermentation time 
(Fig.3.1, 3.2 &3. 3). Again, velocity streamlines has been also solving from naiver stokes equation to see flow condition 


with a time of fermentation. 








Figure 3.1 Temperature profile (inner), and vorticity field at 1.6 m?/hr coolant flow rate 





@IJTSRD | Unique Paper ID-ITJTSRD38009 _ | 


Volume-5|Issue-1 | 


November-December 2020 Page 810 


International Journal of Trend in Scientific Research and Development (IJTSRD) @ www.ijtsrd.com eISSN: 2456-6470 


The results also included vorticity field with fermentation time from naiver stokes equation to see the formation of the 
vortex at a different coolant flow rate. During beer fermentation, the temperature of fermentation is controlled by liquid 
circulated ammonia (Anderson, 2006). Vortex formation during fermentation controls mass transfer and the rate of 
dissolved oxygen (Rubenberger, 2006). This liquid ammonia is controlled by its flow rate. The model was tried to relate 
the coolant flow rate to temperature profile from an energy balance (eqn land 2 above). 


The modeled results were indicating that all phenomena of beer fermentation like air flow (oxygen consumption), 
temperature profile with time, padic pressure with fermentation time and formation of the vortex have influenced by the 
coolant flow rate (Fig.3.1, 3.2 and 3.3). The results indicated that the increasing coolant flow rate decrease temperature 
profile within the fermentation tank (Fig.3.1, 3.2 and 3.3).This indicates that the increasing coolant flow rate increases the 
heat exchange area, hence the high amount of heat absorbed or removed from the coolant from the tank which produced 
exothermically as sugar decomposes to ethanol and carbon dioxide with fermentation time. 


> 
c >i) 3 Trme=21600 s i 


Surface: Temperatureidea 


Temperature (Gee C) 


cS 





Figure 3.2 Temperature profile (inner), and vorticity field at 1.3m?/hr coolant flow rate. 


The velocity streamline was indicating that the flow condition of the fermentation was from bottom to top (upside down) 
at all coolant flow rates and also possible to judge the uniformity of flow at these coolant flowrate, which is better at 1.2 
m3/hr and slightly uniformity decreased for the rest offlow rate (Fig.3.1, 3.2 and 3.3). This indicated that as the coolant 
flow rate increased, the temperature profile within the fermentation tank cannot form uniform mixing which is not 
favorable for mass transfer rate. 


The modeled vorticity results were also indicating that at 1.2 m?/hr coolant flowrate, uniform flow of fluid was observed 
around the stirring center than other coolant rates, which is responsible for a high transfer rate of dissolved oxygen 
(Fig.3.1, 3.2 and 3.3). 


+ . ‘T : 
Ec({l)=1.2 Time=21600 ; 4 Fe({1)=1.2 Time=21600 s 
Surface: Temperature deg ~ . Sstreamiine VOTE) neig 


U2 


Temperature(cegc) 


is 

1é 

7 14 
m S 


- Ud 





Figure 3.3 Temperature profile (inner), and vorticity field at 1.2 coolant flow rate 


3.2. Effect of temperature on fermentation time and beer flavor quality 
The biochemical reaction has been simulating at a different temperature, which already modeled at a different coolant 
flow rate by the finite element method. 


3.2.1. Effect of temperature on nutrient consumption rate. 
Isoluicine was depleted after 130 hours at 12.5°C, 80 hrs at 18°C and 16 hrs at 20°C from the start of fermentation time 
(Fig 3.4, 3.5 and 3.6). The time of consumption at 18 °C and 20°C was decreased by 10 hrs and 30 hrs than at 16 °C (current 
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company value) temperatures respectively. Valine was keeping, decreasing until almost 195 hrs at 12.5°C, 120 hrs at 18°C 
and 85 hrs at 20°C (Fig 3.4, 3.5 and 3.6). 


The amount of Valine in the initial concentration of wort must Keep to low concentration (J. R. Dickinson et al., 2003) to 
maintain higher alcohol production to below thresh hold value. Since the growth of yeast was still maintained after 
depletion of Valine to zero concentration, Valine concentration was counting negative concentration until both maltose 
and malt trios became a zero concentration (Fig 3.4, 3.5 and Fig 3.6). This indicated that the initial concentration of Valine 
is not enough to consume all maltose and malt trios. 


The simulated results of the nutrient were again indicating that the given initial nutrient was become zero after 210 hrs at 
12.5°C, 110 hrs at 18°C and 90 hrs at 20°C for luicine consumption (Fig.3.4, 3.5 and 3.6). These times of consumption were 
indicating that decreasing fermentation temperature increase consumption time of leucine as well as isoleucine. 


3.2.2. Effect of temperature on rate of esters produced. 

The simulated result of esters was indicating that the amount of ethyl acetate formed by fermentation time is relatively 
higher than ethyl caproate and isoamy] acetate in mole per cubic meter (Fig. 3.7, 3.8 and 3.9). Higher alcohols and esters 
are desirable volatile beer constituents, with a few exceptions (Hanke, 2010). Only isoamyl acetate (banana-like aroma) 
concentration is usually above the threshold level in most lager beers, ales normally have ethyl acetate (solvent-like 
aroma) and ethyl hexanoate (apple-like aroma) as supplementary flavoring compounds with levels above their threshold 
taste. Ethyl acetate is the most common ester produced by yeast (Hanke, 2010). 


Several by-products of amino acid metabolism in yeast yield (produce) different flavors in the beer. These flavor-active 
esters are formed by the condensation reaction between either acetyl/Acyl-CoA and higher alcohols or 
ethanol (Warnasoriya, 2010). The final concentrations of higher alcohols and ethyl ester or acetate ester derivatives are 
therefore dependent on the uptake efficiency of the corresponding amino acid and the sugar utilization rate. The simulated 
results were also indicating that ethyl acetate formation follows consumption of Valine amino acids and maltose sugar 
(Fig.3.6, 3.7&3.8). 


From the simulated result, ethyl acetate was reachingit’s thresh hold value of 0.114 Mol/m after 70 hrs at 12.5°C hours, 
30 hrs at 18°C and 22 hrs at 20°C. But, isoamyl acetate maintained thresh hold value of 0.0105 Mol/m (Engan, 2008) until 
the initial concentration of sugar and amino acid were depleted (throughout fermentation time). This indicated that unless 
ethyl ester, isoamyl ester has no flavor effect on the specified condition. 


The simulated results were indicating that esters-based flavor active compounds were obtained at an early stage of 
fermentation time at the specified condition which agrees with actual beer fermentation (Figadu, 2012).Atthe 
fermentation time increased the production rate of ethyl acetate and ethyl caproate was increasing at a decreasing rate. 
This decreasing rate has come from depletion of the initial concentration of both glucose and nutrient. 


Nutrient model concentration [mol/m?] Nutrient model concentration [mol/m*] 


—oO— |soluicine concentration 
—— |uicine concentration 
—4— valine concentration 


—o— |soluicine concentration 
—+— |uicine concentration 
—4— valine concentration 


Concentration [mol/m a 
Concentration [mol/m*} 





0 50 100 150 200 250 0 50 100 150 200 250 


Time [hours] Time [hours] 
Figure 3.4Simulated result of nutrient model Figure 3.5 Simulated result of nutrient model 
concentration at the temperature 12.5 °C, 50 concentration at temperature 18°C, 50 Mol/m? 
Mol/m# pitching rate and 1.3 Mol/m4, 0.6 pitching rate and 1.3 Mol/m4, 0.6 Mol/m? and 2.1 
Mol/m#? and 2.1 Mol/m? concentration (initial) ©Mol/m? concentration for isoluicine, luicine and Valine 
for luicine, isoluicine and Valine respectively. respectively. 
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Nutrient model concentration [mol/m*] 


—0— |soluicine concentration 
—— |uicine concentration 
—4— valine concentration 


Figure 3.6Simulated result of nutrient model 
concentration at temperature 20°C, 50 Mol/m? 
pitching rate and 1.3 Mol/m4, 0.6 Mol/mé? and 2.1 
Mol/m# concentration (initial) for isoluicine, luicine 
and Valine respectively. 


Concentration [m ol/m>] 





0 50 100 150 200 250 
Time [hours] 


Ethyl acetate concentration [mol/m*] 


Ethyl acetate concentration 


Esters formation concentration [mol/m?*] 


0.014 
—o— Ethyl caproate concentration 
—+— |soamy! acetate concentration | 
0.012 = 
0.01 


0.008 


Concentration [mol/m*} 


0.002 





oa = 
0 50 100 150 200 250 0 50 100 150 200 250 


Time [hours] Time [hours] 


Figure 3.7 Esters model simulated result at temperature 12.5 °, pitching rate 50 Mol/m? and 200 
Mol/m?, 320 Mol/m? and 100 Mol/m concentration (initial) for glucose, maltose and malt triose 


respectively. 
Ethyl acetate concentration [mol/m?] Esters formation concentration [mol/m*] 
0.018 
—°©— Ethyl caproate concentration 
0.016 ——* = Isoamyl acetate concentration 
0.014 7 
“o 
£ 0.012 
© 2 
E = 0.01 
0 3 
é E 
< = 0.008 
Ww S 
e 
6 0.006 
oO 


0.004 


0.002 





0 50 100 150 200 250 0 50 100 150 200 250 
Time [hours] Time [hours] 


Figure 3.8 Esters model simulated result at temperature 18 °, pitching rate 50 Mol/m? and 200 Mol/m?, 
320 Mol/m? and 100 Mol/m? concentration (initial) for glucose, maltose and malt triose respectively. 
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: 3 
Ethyl acetate concentration [mol/m*] Esters formation concentration [mol/m*] 


0.7 0.018 





—o2— Ethyl caproate concentration 
06 0.016 —+— |soamy! acetate concentration 

0.014 
0.5 a 

£ 0.012 
“a 2 
0.4 = 

E — 0.01 
= 2 
= 0 

< — 0.008 
< 0.3 © 
eS 

oO 0.006 
0.2 2 

0.004 

0.1 0.002 

0 oe 
0 50 100 150 200 250 0 50 100 | 150 200 250 
Time [hours] Time [hours] 


Figure 3.9Esters model simulated result at temperature 20 °C, pitching rate 50 Mol/m? and 200 Mol/m3, 
320 Mol/m? and 100 Mol/m? concentration (initial) for glucose, maltose and malt triose respectively. 


Biomass concentration [mol/m 4 


yeast concentration 


Biomass concentration [mol/m ") 


yeast concentration 


160 


140 


=_k 
hm 
Oo 


X [mol/m?] 
X [mol/m?] 


60 





40 


0 50 100 150 200 250 ? on _ a —_ = 
Time {hours} Time [hours] 
Figure 3.10Simulated result of yeast at 12.5°C, Figure 3.11Simulated result of yeast and ethanol at 
50 Mol/m? pitching rate and 1.3 Mol/m4, 0.6 18°C, 50 Mol/m*pitching rate and 1.3 Mol/m3, 0.6 
Mol/m? and 2.1 Mol/m#? concentration for Mol/m3 and 2.1 Mol/m3 concentration for isoluicine, 
isoluicine, luicine and Valine respectively. luicine and Valine respectively. 


Biomass concentration [mol/m 34 


east comce nitration 


0 Oh 





150 


X {malin 


oo 


=O 
Ch =O} 100 150 CO eC} 


Tinne [hours] 


Figure 3.12Simulated result of yeast and ethanol at 20°C, 50 Mol/mé pitching rate and 1.3 Mol/m3, 0.6 
Mol/m3 and 2.1 Mol/m3 concentration for isoluicine, luicine and Valine respectively. 


4. MODEL VALIDATION 

To see the accuracy of the generated data, validation of the model had been necessary by comparing with already existed 
data or experimental data. In this case, the experiment was done to compare experimental value with a numerical value. 
Accordingly, for heat transfer model, temperature profile was collected along the arc length of the tank of the experimental 
Staff at 1.4 coolant flow rate (current coolant flow rate) and compared with the modeled results as follows. 
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nparison of simulated and reference data at 16 degC, 1.4 coolant flowrate and 21600s with arc length 3 


Temperature(degt) 


0 0,05 0.1 0.15 0.2 
Arc length(m) 


Figure 4.1 Comparison of simulated and experimental data at 16°C, 1.4m?2/hr coolant flow rate and 21600 second 
along arc length 


5. Conclusion 

The beer fermentation process has been modeled by 
having a phenomenon that can be representative of beer 
fermentation. Energy balance was done to have a 
temperature profile with fermentation time. Relevant 
kinetics were developed for yeast growth, nutrients and 
esters model and simulated by incorporating operating 
parameter like temperature. 


The finite element method was used to solve energy 
balance that has solved from an energy balance at a 
different coolant flow rate of ammonia coolant. The 
modeled result was indicating that as the coolant flow rate 
decreased, the temperature profile with fermentation time 
was increased. 1.2 m3/hr, 1.3 m?/hr, and 1.6 m?/hr were 
the selected coolant flow rates to see temperature profiles 
with fermentation time. This is due to of that, most of the 
brewery share company used 1.4 m%/hr coolant flow rate 
currently which the study wants to vary up and down 
from the current one. The temperatures that were 
observed at 1.2 m?/hr, 1.3 m?/hr and 1.6 m?/hr coolant 
flow rate were 20°C, 18°C and 12.5°C respectively with 
fermentation time. 


The modeled kinetics of flavors (esters), nutrients and 
growth (yeast) model were simulated with the above 
temperatures to see its effect on fermentation time and 
beer flavor quality. The simulated results were indicated 
that in 20°C time to reach a flavor threshold value and 
time of Valine, luicine and isoluicine consumption 
wereshorter than the other two temperatures. 
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LIST OF ABREVIATION 
E Ethanol concentration (mol/m?) 
FE Finite Element 
G Glucose concentration (mol/m?) 
IA Isoamyl alcohol concentration (mol/m?) 
IAc Isoamyl acetate concentration (mol/m?) 
IB Isobutyl alcohol concentration (mol/m+?) 
M Maltose concentration (mol/m?) 
MATLAB Matrix Laboratory 
MB 2-methyl-1-butanolconcentration (mol/m+?) 
N Malt triose concentration (mol/m+?) 
P n- propanol concentration (mol/m+?) 
LIST OF SYMBOLS 
AH¢c Over all heat of formation for glucose (kJ/mol) 
Li Specific rate of glucose uptake (h‘!) 
[2 Specific rate of maltose uptake (h‘') 
[3 Specific rate of malt triose uptake (h*!) 
Ui Maximum reaction velocity for i sugar (i= G, M or N, h*?) 
Ly Specific rate of isoluicine uptake (h*!) 
Li Specific rate of luicine uptake (h‘‘) 
La Specific rate of isoamyl alcohol formation (h‘!) 
Lio Arrhenius frequency factor for maximum velocity (h*!) 
[Lx Specific rate of yeast growth (h‘!) 
A Area of heat transfer (m2) 
Cp Specific heat capacity of wort (kJ/kg °C) 
D First order time delay (Michiko, Hiroshi, & Suteaki) 
Eui Arrhenius activation energy for maximum velocity (cal/mol) 
E'pi Arrhenius activation energy for inhibition constant (cal/mol) 
Eig Arrhenius activation energy for michaelis constant (cal/mol) 
Kaal Effective first order rate constant for uptake of AAI (m?/mol h) 
Ki Michaelis constant for it* sugar (i = G, M or N, mol/m*) 
K'i Michaelis constant for it sugar (i = G, M or N, mol/m?) 
Kio Arrhenius frequency factor for michaelis constant (mol/m/3) 
K'io Arrhenius frequency factor for inhibition constant (mol/m?) 
Kx Empirical yeast growth inhibition constant (mol/m+?) 
Q Rate of heat flow (kJ/mol h) 
R Gas constant (cal/mol K) 
t Time(Michiko et al.) 
T Temperature (°C) 
U Heat transfer coefficient (W/m2°C) 
X Biomass concentration (mol/m+?) 
Xo Initial yeast concentration (mol/m?) 
YEG Yield coefficient, mole E per Mole G 
YEM Yield coefficient, mole E per Mole M 
YEN Yield coefficient, mole E per Mole N 
Yxc Yield coefficient, mole X per Mole G 
Yxm Yield coefficient, mole X per Mole M 
YXn Yield coefficient, mole X per Mole N 
0 Density of wort (kg/m) 
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